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Abstract 

We present the Ernst potential and the Hne element of an exact solution of Einstein's vacuum 
field equations that contains as arbitrary parameters the total mass, the angular momentum, and 
the quadrupole moment of a rotating mass distribution. We show that in the limiting case of slowly 
rotating and slightly deformed configuration, there exists a coordinate transformation that relates 
the exact solution with the approximate Hartle solution. It is shown that this approximate solution 
can be smoothly matched with an interior perfect fluid solution with physically reasonable prop- 
erties. This opens the possibility of considering the quadrupole moment as an additional physical 
degree of freedom that could be used to search for a realistic exact solution, representing both the 
interior and exterior gravitational field generated by a self-gravitating axisymmetric distribution 
of mass of perfect fluid in stationary rotation. 
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I. INTRODUCTION 



Astrophysical compact objects are in general not spherically symmetric rotating mass dis- 
tributions. To describe the corresponding gravitational field one can assume axial symmetry, 
with an axis of symmetry that coincides with the axis of rotation. Then, the deviation from 
spherical symmetry can be described by means of axisymmetric multipole moments. 

In general relativity, vacuum solutions with multipole moments have been known for 
a long time. In fact, if we restrict ourselves to solutions with only monopole moment, 
Birkoff's theorem guarantees that the only solution is the Schwarzschild metric. Moreover, 
if rotation is also taken into account, the black hole uniqueness theorems state that the only 
asymptotically fiat solution with a regular horizon is the Kerr metric The next interesting 
multipole is the quadrupole. In this case, the uniqueness theorems do not apply and it is 
possible to find a large number of different vacuum solutions with the same quadrupole. 
Differences appear only at the level of higher multipoles. The first static solution with 
an arbitrary quadrupole was found by Weyl using cylindrical coordinates. Later on, 
Erez and Rosen (ER) discovered a solution with arbitrary quadrupole in prolate spheroidal 
coordinates which are more convenient for the investigation of multipole solutions. Zipoy 
3 1 and Voorhees j4| found a simple transformation which allows to generate static solutions 
from a given one. In particular, the Zipoy- Voorhees (ZV) transformation can be used to 
generate the simplest solution with quadrupole, starting from the Schwarzschild metric. 
Stationary solutions represent an additional chal 
rotating solution was found only in 1963 by Kerr 
axisymmetric fields was an important achievement that allowed to search for the symmetries 
of the field equations upon which modern solution generating techniques are based. 



enge. In fact, the first physically relevant 
^. The Ernst representation of stationary 



In this work, we will limit ourselves to the study of a particular stationary axisymmetric 
vacuum solution which was derived by Quevedo and Mashhoon (QM) in [?, ^ as a gener- 



alization of the ER metric [2]. The QM solution contains in general an infinite number of 
gravitational and electromagnetic multipole moments. Here, however, we neglect the elec- 
tromagnetic field and focus on the contribution of the gravitational quadrupole only. The 
main goal of the present work is to show that the QM solution can be used to describe the 
exterior gravitational field of rotating compact objects. First, we will see that the set of 
independent and arbitrary parameters entering the metric determines the total mass, an- 
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gular momentum and mass quadrupole moment of the source. Although higher multipole 
moments are present, they all can be expressed in terms of the independent lower multipoles. 
This result is obtained by using the invariant definition of relativistic multipole moments 
proposed by Geroch and Hansen ^, [l^ . The spacetime turns out to be asymptotically flat 
and and free of singularities outside a region which can be "covered" by an interior perfect 
fluid solution. This last property, however, is shown only in the limiting case of a slightly 
deformed body with uniform and slow rotation. 

There exists in the literature a reasonable number of interior spherically symmetric solu- 
tions which can be matched with the exterior Schwarzschild metric. Nevertheless, a major 
problem of classical general relativity consists in finding a physically reasonable interior so- 
lution for the exterior Kerr metric. Although it is possible to match numerically the Kerr 
solution with the interior field of an infinitely tiny rotating disk of dust 17|, such a hypo- 



thetical system does not seem to be of relevance to describe astrophysical compact objects. 
It is now widely believed that the Kerr solution is not appropriate to describe the exterior 
field of rapidly rotating compact objects. Indeed, the Kerr metric takes into account the 
total mass and the angular momentum of the body. However, the moment of inertia is an 
additional characteristic of any realistic body which should be considered in order to cor- 
rectly describe the gravitational field. As a consequence, the multipole moments of the field 
created by a rapidly rotating compact object are different from the multipole moments of 
the Kerr metric. For this reason a solution with arbitrary sets of multipole moments, such 
as the QM solution, can be used to describe the exterior field of arbitrarily rotating mass 
distributions. 

To completely characterize the spacetime it is necessary to find an exact solution with 
a set of interior multipole moments. Due to the generality of the exterior solution, one 
can expect that the corresponding interior solution can be derived by postulating an infinite 
series of inner multipoles with arbitrary functions which can be matched one by one with the 
exterior solution. To see if this procedure is realizable, in this work we analyze the special 
case of a deformed rotating body with quadrupole moment only. Despite this simplification, 
the problem of finding interior solutions is at present still out of reach due, in part, to the 
complexity of Einstein's equations with a realistic model for the inner configuration. We 
therefore limit ourselves here to the study of approximate perfect fluid solutions. 

In fact, in the case of slowly rotating compact objects it is possible to find approximate 
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interior solutions with physically meaningful energy-momentum tensors and state equations. 
Because of its physical importance, in this work we will study the Hartle-Thorne (HT) 
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13| interior solution which can be coupled to an approximate exterior metric. Hereafter 
this solution will be denoted as the HT solution. One of the most important characteristics 
of this family of solutions is that the corresponding equation of state has been constructed 
using realistic models for the internal structure of relativistic stars. Semi-analytical and 
numerical generalizations of the HT metrics with more sophisticated equations of state have 
)een proposed by different authors. A comprehensive review of these solutions is given in 



In all these cases, however, it is assumed that the multipole moments (quadrupole 
and octupole) are relatively small and that the rotation is slow. We will find the explicit 
coordinate transformation that transforms an approximation of the QM solution into the 
exterior HT solution and can be matched with an approximate solution. In this manner, 
we show that the QM solution satisfies the main physical conditions to be matched with an 
interior solution. 

II. EXTERIOR SOLUTION 



It is well known that the main physical information about an axisymmetric stationary 
vacuum gravitational field can be encoded in the complex Ernst potential jo]. For a special 
case of the QM solution the Ernst potential reads 

^ ^ x-l-{x + 1)X^ + i[y{X + ^) + X- fi] 

where the function ijj can be written in terms of the Legendre polynomials P and Legendre 
functions Q as 

^ = (1 - 6)Qo - 6qP2Q2 , P2 = liSy' - 1) , (2) 

^ ^ X. X + 1 ^ 1,„2 1NT^ + 1 3 

Qo = - in , CJi = — m x , U2 = -[ox — 1) m x . (6) 

Here 6 and q are arbitrary real constants. Moreover, A and /i are functions of x and y 

A = a{x' - lY-\x + yf'-^e'''^^- , ^ = a{x' - l^'^x + yf^-'e'^'^^^ , (4) 



with 



lln^-^^ + l{l-y'Txy) + l[x{l-y^)Ty{x^-mn^^ 



The constant a can be represented in terms of the additional constants a and m as 

a — m 



a 



a 



a = V m? — . 



(6) 



The solution is asymptotically flat and possesses the following independent parameters: m, a, 
5, and q. In the limiting case a = 0, a = 0, g = and 5 = 1, the only independent parameter 
is m and the Ernst potential ([T]) determines the Schwarzschild spacetime. Moreover, for 
a = a = and g = we obtain the Ernst potential of the ZV static solution which 
is characterized by the parameters m and 5. Furthermore, for a = a = and 5 = 1, 
the resulting solution coincides with the ER static spacetime [2)]. The Kerr metric is also 
contained special case for g = and 5 = 1. 

The physical significance of the parameters entering the potential ([1]) can be established 
in an invariant manner by calculating the relativistic Geroch-Hansen multipole moments. 

n 

We use here the procedure formulated in pj| which allows us to derive the gravitoelectric 
Mn as well as the gravitomagnetic Jn multipole moments. A lengthly but straightforward 
calculation yields 

M2fc+i = J2fc = 0, fc = 0,1,2,... 
Mo = m + (t(5 - 1) 
35^ - 45 + 6) - ma'^5{5 - 2) - 3mV(5 - 1) - 
Ji = ma + 2acr(5 — 1) , 



M2 = —a^ 
15 



3 ^ 



J3 = — aa'^5q — a 
^ 15 ^ 



3 ^ 



352 - 5 + 3) + ma\35^ - 65 + 2) + 4mV(5 - 1) 



m 



(7) 
(8) 
(9) 
(10) 

(11) 



The even gravitomagnetic and the odd gravitoelectric multipoles vanish identically because 
the solution possesses and additional reflection symmetry with respect to the hyperplane 
y = 0. Higher odd gravitomagnetic and even gravitoelectric multipoles can be shown to be 
linearly dependent since they are completely determined in terms of the parameters m, a, q 
and 5. From the above expressions we see that the ZV j^, 4| parameter changes the value of 
the total mass Mq as well as the angular momentum Ji of the source. The mass quadrupole 
M2 can be interpreted as a nonlinear superposition of the quadrupoles corresponding to 
the ZV, ER and Kerr spacetimes. Notice that in the limiting static case of the ZV metric 
the only non-vanishing parameters are m = a and 5 so that all gravitomagnetic multipoles 
vanish and we obtain Mq = m5 and M2 = 6m'^{l — 5^)/3 for the leading gravitoelectric 
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multipoles. This means that the ZV solution represents the field of a static deformed body. 
To our knowledge, this is the simplest generahzation of the Schwarschild solution which 
includes a quadrupole moment. 

In order to completely describe the geometric and physical properties of the spacetime, 
it is convenient to calculate the explicit form of the metric. In fact, the Ernst potential is 
defined as 



f + in , with a{x^ - l)n^ = fcuy , a{l - y^)ny = -f 



(12) 



where / and co are the main metric functions which determine the line element in prolate 
spheroidal coordinates {t,x,y,ip): 



ds^ = fidt-ud^y - 



J 



dx^ 



(a;2-l)(l-l/W 



(13) 



^x'^ — 1 1 — y"^ J 

The metric function 7 can be calculated by quadratures once / and ou are known. The 
calculation of the metric functions results in 



/ 



R 



-2qSP2Q2 



M 



^ ^ _2a - 2(7^e2«'^^^«^ 
R 



where 



a J {x^-iy 

R = a+a_ + L = + 6^, 

M^(x + ly-^ \x{l - y^)(X + r])a+ + y(x^ - 1)(1 - Xv)b^ 



(14) 

(15) 
(16) 

(17) 
(18) 



1,_ ,2, ^'-1 



7 = -(l + g)^ln' 



2 ^ ' x^ — y^ 



+ 2q{l - P,)Q^ + q\l - + P^ml - Ql) (19) 



+-(x' - l)(2g^ _ ?,xQiQ2 + ?>QoQ2 - Q'2)\- 



(20) 



Furthermore 



a± = (x ± \)'~\x{\ - Xri) ± (1 + Ar^)], (21) 

h^ = {x±\f-\y{\-rri)^{\-ri)\ . (22) 

It is easy to show that this solution is asymptotically fiat and free of singularities outside 
the sphere x — \ which represents a naked singularity. In fact, a numerical investigation 
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FIG. 1: Structure of the hypersurface x = 1 in the spacetime of a rotating mass with arbitrary 
quadrupole parameter. The plot (a) illustrates the case of vanishing quadrupole, q = 0, with 
6 = 1. The hypersurface x = 1 corresponds to an event horizon (dashed line). The horizon covers 
the ring singularity which is caused by the rotation of the body. The case (b) is for 5 = 2 and 
q = 10. The horizon x = 1 becomes singular due to the presence of the quadrupole moment so 
that the spacetime is characterized by a naked singularity (solid line). By varying the value of 
the quadrupole parameter q it is possible to generate additional singularities (c) which are always 
inside the outer naked singularity situated at x = 1. 

of the Kretschmann invariant shows that the hypersurface x = 1 is singular, independently 
of the value of the parameters m, 6 1, a, and q. Only in the limiting case g = and 
6 = 1, the hypersurface a; = 1 coincides with the exterior horizon of the Kerr spacetime. 
This situation is illustrated in Fig{T] which shows that a naked singularity is always present 
when the quadrupole parameter q is non zero or when 6 1. Moreover, the symmetry axis 
y = ±1 is free of singularities (except at x = 1), indicating that the condition of elementary 
flatness is satisfied 

All the properties mentioned above seem to indicate that the solution can be used to 
describe the exterior field of a deformed, rotating mass distribution. From the point of view 
of general relativity, however, to describe the entire manifold of a rotating body it is necessary 
to find an interior solution which takes into account the inner properties of the body and 
that can be matched to the exterior vacuum solution. However, rather few exact stationary 
solutions that involve a matter distribution in rotation are to be found in the literature. In 
particular, the interior solution for the rotating Kerr solution is still unknown. In fact, the 
quest for a realistic exact solution, representing both the interior and exterior gravitational 
field generated by a self-gravitating axisymmetric distribution of mass of perfect fiuid in 
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stationary rotation is considered a major problem in general relativity. We believe that the 
inclusion of a quadrupole in the exterior and in the interior solutions adds a new physical 
degree of freedom that could be used to search for realistic interior solutions. To see if this 
is true, we will use a special limit of the above solution with quadrupole moment and will 
show that in fact it can be matched with a realistic approximate interior solution. 



III. THE FIELD OF A SLOWLY ROTATING AND SLIGHTLY DEFORMED 
BODY 



In this section we calculate the limiting case where the deviation from spherical symmetry 
is small and the body is slowly rotating. It is convenient to introduce the new spatial 
coordinates r and 6 by means of {igI 



X 



r — m 
a 



y 



cos 9 



(23) 



and to choose the ZV parameter as 5 = 1 + sg, where s is a real constant. Then, expanding 
the metric ( |T3l) to first order in the quadrupole parameter q and to second order in the 
rotation parameter a, we obtain 



/ =1 
+3g 



2m 2a^mcos^^ , , / 2m , , 
+ ^ + g(l + s) 1 ) In 



r 
r 



^ 2m\ 
r / 



1 



1 



m 
r 



3 cos^ 9 - 1 



2m 

+ ^ f _!L _ 1^ (3 cos^ 0-l)-j sin^ In (^1 



\2m 



2m \ 
r / 



2amr sin^ 9 
r — 2m ' 



(24) 
(25) 



1 , r(r — 2m) 

7 = — In 

2 (r — mY — To^ cos^ 9 



rn? cos^ 9 sin^ 9 



r{r — 2m) ((r — m)^ — m? cos^ i 



+q{l + s) In 



r(r — 2m) 



(r — m)2 — im? cos^ 9 



-3q 



1 / r \ , / 2m 
1 + - l lnl 

2 \m J \ r 



sin^ 9. 



(26) 



Furthermore, we introduce coordinates R and O by means of 



R + Mq + -A^gsin^ 6 



R__ 

M ^ 2M^ 



1 - 



2M 



1 



2M 



R 



a 

2R 



1 + 



2M 



R 



M\ 27W\ / 3M\ 



(27) 



6' = 6 — sin cos < -o 
1 2^ 



2 + 



R 
M 



1 In 1 



2M\ 



where 



-ma, Q = m q . 

m 5 



M = m(l - q), J 

A straightforward computation shows that the metric (!24l) - (!26l) can be written as 
2M\ 

' I ^ 'V R J R" 



(28) 
(29) 



R J 



1 



2M\ 
If) 



l-2A;2P2(cos0)-2 1- 



R 
J 



R^ 



dR' 



-R^[l - 2k3P2{cos 0)](d02 + sin^ 0#^) + 4- sin^ Qdtd(j) 



R 



(30) 



where 



k2 
k3 



J2 



MR^ V i? ; 8 ^'\M 



ki- 



6J2 
i?4 



j2 _5g-JV7W / 



27W 



-1/2 . D 



and we choose for the free parameter s the particular value s = — 1 for the sake of simplicity. 
Here QY^ are the associated Legendre functions of the second kind 



Qlix) 



1) 



1/2 



3x^-2 



3 , x+ 1 

-xm 

2 x-1 



, Qlix) = -{x^-l)\n 



x + 1^5x-3x3^ (31) 
X — 1 x^ — I 



The line element f pOj) was first derived in an equivalent representation by Hartle and Thorne 
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13[. Consequently, our results show that the general solution (fT3l) -(j22l) contains as a 



special case the HT metric which describes the exterior gravitational field of any slowly and 
rigidly rotating deformed body. 

In the case of ordinary stars, such as the Sun, the smallness of the parameters 



un 



lo-^ Js.n/7^; 



metric f l30l) to obtain 

ds" - 



18| 



Sun 



10 , Qsun/T^Sun 



dt^ + ^ sin^ Qdtdcj) 



10 allows us to simplify the 
2M 2Q 



1 



2Q 



P2(COS0) 



dR^ 
(32) 



The accuracy of this metric is of one part in 10^^. Consequently, it describes the gravitational 
field for a wide range of compact objects, and only in the case of very dense {Ai ~ TZ) or 
very rapidly rapidly rotating (J ~ 71'^) objects large discrepancies will appear. 



IV. THE INTERIOR SOLUTION 



If a compact object is rotating slowly, the calculation of its equilibrium properties reduces 
drastically because it can be considered as a linear perturbation of an already-known non- 
rotating configuration. A very good approximation of the internal structure of the body is 
delivered by a perfect fluid model satisfying a one-parameter equation of state, V = V{£), 
where V is the pressure and £ is the density of total mass-energy. A further simplification 
follows from the assumption that the configuration is symmetric with respect to an arbitrary 
axis which can be taken as the rotation axis. Furthermore, the rotating object should be 
invariant with respect to reflections about a plane perpendicular to the axis of rotation, 
i.e, about the equatorial plane. As for the rotation, it turns out that configurations which 
minimize the total mass-energy (e.g., all stable configurations) must rotate uniformly 14 1. 
Uniform rotation facilitates the study of the internal properties of the body, especially if we 
restrict ourselves to the case of slow rotation. That is, we assume that angular velocities Vl 
are small enough so that the fractional changes in pressure, energy density and gravitational 
field due to the rotation are all less than unity, i.e. Vt^ <^ (^^^ where Ai is the mass and 
TZ is the radius of the non-rotating configuration. The above condition is equivalent to the 
physical requirement Q <^ c/TZ. The metric functions must be found by solving Einstein's 
equations 

R; - l^^^R = Svtt; , (33) 

where the stress-energy tensor is that of a perfect fluid 

t; = {s + v)u^u, - V5;. (34) 

The 4-velocity which satisfies the normalization condition u^u^ = 1 is 

«« = M0 = O, u^ = ^lu\ u' ={gu + 2ngt^ + n''g^^y^'\ (35) 

where the angular velocity is a constant throughout the fluid. 

The explicit integration of the inner Einstein's equations depends on the particular choice 
of the density function S{R), where R is the radial coordinate, and on the equation of state 
V = V{£). For the sake of simplicity, we consider here only the simplest case with S =const 
and let the pressure V be determined by the field equations. Then, the total mass of the 
non-rotating body is = inSTl^/S, where 71 is the radius of the body. Under these 
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conditions, the resulting line element can be written as 

MR) 



[I + 2$) dr - 



l + 2i?- 



dR 



+ $2(i?)P2(cose) 



dR' 



where 



-R^ [1 + 2$2(i?)i^2(cos e)] [dO^ + sin^ 0(^0 - udtf], 



$ = $o(i?) + $2(i?)P2(cose), 



(36) 



(37) 



is the interior Newtonian potential. Here $o is the interior Newtonian potential for the 
non-rotating configuration and ^2{R) is the perturbation due to the rotation. Moreover, 



H) = VL — ui 



(38) 



is the angular velocity of the local inertial frame, where =const is the angular velocity of 
the fluid and uj — Cu{R) is the angular velocity of the fluid relative to the inertial frame. With 
a suitable choice of the integration constant, the interior unperturbed Newtonian potential 
can be expressed as 

MR) = -27rp - ^] , (39) 



where TZ is the radius of the non- rotating conflguration. The function ^2{R) satisfles the 
following equation 



d^R) _ ^ (r.. (^^R'P 2 \ 2x 47r , , 



where x is deflned by 



_ 2M Stt 2o3 

dR - -^^'^^^ + yP^^ ' 



(40) 



(41) 



For the angular velocity relative to the local inertial frame u one obtains the differential 

equation 



1 d 



. du 



4 . u,Lv I 4 dj 

^ZrI ■'m^RdR^ 



0, 3{R) 



4vr „n 



1/2 



(42) 



which is related to the equilibrium condition of a rotating self-gravitating body, namely, 
the condition that there exists a balance between pressure forces, gravitational forces and 
centrifugal forces. The solution for u{R) must be regular at the origin and outside the body 
it must take the form (jj[R) = Vt — 2J jR?^ where J is the total angular momentum of the 
star. This condition guarantees a smooth matching of the corresponding metric functions 
on the surface of the mass distribution. 
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For the integration of the field equations it is convenient to consider also an additional 
differential equation which follows from the conservation law T^'^^^ = and involves the 
pressure and its partial derivatives. Even in the simple case of £ = const, the resulting 
expression is rather cumbersome. A detailed analysis of this equation will be presented 
elsewhere 18|]. The integration of the system of partial differential equations (HU]) - 



together with the differential equation for the pressure, cannot be carried out analytically. 
Nevertheless, it is possible to find numerical solutions not only in the simple case S = const, 
but also in the case of more realistic density functions jisl . 

It is easy to show analytically that the interior solution f l5B]) can be matched smoothly 
on the surface R = TZ with the exterior solution f l30p in the special case of an unperturbed 
configuration with $2 = and $o as given in Eq. fl59]) . In the general case of a rotating 
deformed body, the matching can be performed only numerically. In fact, the solutions 
of the field equations for the inner distribution of mass are calculated using the matching 
conditions as boundary conditions. In this manner, it is possible to find realistic solutions 
that describe the interior and exterior gravitational field of slowly rotating and slightly 
deformed mass distributions. 



V. CONCLUSIONS 



In this work we studied an exact solution of Einstein's vacuum field equations which 
represents the exterior gravitational field of a stationary axisymmetric rotating mass distri- 
bution. It contains four free parameters that determine the total mass, angular momentum, 
and quadrupole momentum of the body. It was shown that the solution is asymptotically 
fiat and is free of singularities outside a region that can be covered by an interior solution. 
In fact, we show that the presence of a gravitoelectric quadrupole changes drastically the 
geometric properties of the spacetime. Independently of the value of the quadrupole, there 
always exists a singular surface which is not surrounded by an event horizon. This implies 
that the gravitoelectric quadrupole moment of the source can always be associated with the 
presence of naked singularities. However, the naked singularities are situated very close to 
the origin of coordinates so that a physically reasonable interior solution could be used to 
cover them. 

The main goal of this work was to take the first step to prove that the exact exterior 
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QM solution can be matched in general with an interior solution. To this end we calculated 
the limiting case of an approximate solution for a slowly rotating and slightly deformed 
source. This means that the solution is expanded to first order in the quadrupole parameter 
and to second order in the rotation parameter. We show that a particular choice of the 
ZV parameter allows us to find explicitly a coordinate transformation that transforms the 
solution into the approximate exterior Hartle metric. Assuming that the inner configura- 
tion can be described by a perfect fluid which rotates uniformly with respect to an axis 
that coincides with the axis of symmetry and that the source is symmetric with respect to 
reflections about the equatorial plane, the set of differential equations which follow from 
Eintein's interior equations was derived and investigated. A particular analytic sohition was 
presented that can be matched with the approximate exterior solution in the case of an un- 
perturbed non-rotating perfect fluid. To consider the general case of a slowly rotating and 
slightly deformed perfect fluid it is necessary to perform numerical calculations to integrate 
the corresponding set of differential equations. Our numerical results show that it is possible 
to obtain numerical solutions with realist physical properties and to match them with the 
exterior approximate solution. The inclusion of the quadrupole parameter in the interior 
and exterior metrics is equivalent to adding a new degree of freedom which facilitates to 
search for interior solutions. 

We conclude that the particular QM solution presented in this work can be used to 
describe the gravitational fleld of a rotating body with arbitrary quadrupole moment. In 
this work, it was shown that in the case of a slightly deformed body with uniform and slow 
rotation, there exist a set of solutions of Einstein solutions which can be used to completely 
describe the spacetime. The case of an arbitrarily rotating perfect fluid with arbitrary 
quadrupole parameter will be investigated in a future work. 
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